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Abstract 

We formulate a dynamical fluctuation theory for stationary non equilibrium states (SNS) which 
is tested explicitly in stochastic models of interacting particles. In our theory a crucial role is 
played by the time reversed dynamics. Within this theory we derive the following results: the 
modification of the Onsager-Machlup theory in the SNS; a general Hamilton-Jacobi equation 
for the macroscopic entropy; a non equilibrium, non linear fluctuation dissipation relation valid 
for a wide class of systems; an H theorem for the entropy. We discuss in detail two models of 
stochastic boundary driven lattice gases: the zero range and the simple exclusion processes. In 
the first model the invariant measure is explicitly known and we verify the predictions of the 
general theory. For the one dimensional simple exclusion process, as recently shown by Derrida, 
Lebowitz, and Speer, it is possible to express the macroscopic entropy in terms of the solution 
of a non linear ordinary differential equation; by using the Hamilton-Jacobi equation, we obtain 
a logically independent derivation of this result. 
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1. Introduction 

The Boltzmann-Einstein theory of equilibrium thermodynamic fluctuations, as described 
for example in Landau-Lifshitz [J20|| , states that the probability for a fluctuation from 



equilibrium in a macroscopic region of volume V is proportional to 

exp{VAS/k} 

where AS* is the variation of entropy density calculated along a reversible transformation 
creating the fluctuation and k is the Boltzmann constant. This theory is well established 
and has received a rigorous mathematical formulation in classical equilibrium statistical 
mechanics via the so called large deviation theory pi ]. The rigorous study of large 



deviations has been extended to hydrodynamic evolutions of stochastic interacting particle 
systems [19]]. In a dynamical setting one may ask new questions, for example what is 
the most probable trajectory followed by the system in the spontaneous emergence of a 
fluctuation or in its relaxation to equilibrium. The Onsager-Machlup theory |24| gives 



the following answer under the assumption of time reversibility. In the situation of a 
linear macroscopic equation, that is, close to equilibrium, the most probable emergence 
and relaxation trajectories are one the time reversal of the other. Developing the methods 
of [[HJ, this theory has been extended to nonlinear hydrodynamic regimes ||17|| . 



In the present paper we formulate a general theory of large deviations for irreversible 
processes, i.e. when detailed balance does not hold. This question was previously ad- 
dressed in []II| for finite dimensional diffusions and in || for lattice gases. Natural exam- 
ples are boundary driven stationary non equilibrium states (SNS), e.g. a thermodynamic 
system in contact with two reservoirs. In such a situation there is a flow of matter or other 
physical property like heat, charge,... through the system. As we shall see, the sponta- 
neous fluctuations of the process are described by the time reversed dynamics, which is 
defined below. 

Spontaneous fluctuations, including Onsager-Machlup symmetry, have been observed 
in stochastically perturbed reversible electronic devices ||22|| . In their work, these authors 
study also non reversible systems and observe violation of Onsager-Machlup symmetry; 
in the present work we shall connect such violations to the time reversed dynamics. 

We are interested in many body systems in the limit of infinitely many degrees of 
freedom. The basic assumptions of our theory are the following. 

1) The microscopic evolution is given by a Markov process X t which represents the 
configuration of the system at time t. This hypothesis probably is not so restrictive 
because also the Hamiltonian case discussed in |7j in the end is reduced to the analysis of 
a Markov process. The stationary non equilibrium state (SNS) is described by a stationary, 
i.e. invariant with respect to time shifts, probability distribution P st over the trajectories 
oiX t . 

2) The system admits also a macroscopic description in terms of density fields which 
are the local thermodynamic variables. For simplicity of notation we assume there is 
only one thermodynamic variable p. The evolution of the field p = p(t, u) where u is the 
macroscopic space coordinate, is given by a diffusion type hydrodynamical equation of 
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the form 

d tP = l -V ■ (D(p)Vp) = \ *h(DM*hP)=*>(p) (1-1) 

l<i,j<d 

The interaction with the reservoirs appears as boundary conditions to be imposed on 
solutions of ( |1 . 1|) . We assume that there exists a unique stationary solution p of i.e. 
a profile p(it), which satisfies the appropriate boundary conditions such that 2)(p) = 0. 
This holds if the diffusion matrix Dij(p) in (|1.1|) is strictly elliptic, namely there exists a 
constant c > such that -D(p) > cl (in matrix sense). 

This equation derives from the underlying microscopic dynamics through an appro- 
priate scaling limit. The hydrodynamic equation ( |1.1| ) represents a law of large numbers 
with respect to the probability measure P s t conditioned on an initial state Xq. The initial 
condition for ( |1 . 1| ) is determined by Xq. Of course many microscopic configurations give 
rise to the same value of p(0, it). In general p = p(t, u) is an appropriate limit of a pN{X t ) 
as the number N of degrees of freedom diverges. 

3) Let us denote by 9 the time inversion operator defined by 9X t = X_ t . The prob- 
ability measure P* t describing the evolution of the time reversed process X£ is given by 
the composition of P s t and 9~ l that is 

p; t {x; = (fH,te [h,t 2 ]) = p st(x t = 4>-ute [-t 2 , -tj) (1.2) 

Let L be the generator of the microscopic dynamics. We remind that L induces the 
evolution of observables (functions on the configuration space) according to the equation 
d t Ex [f(Xt)] = Ex [(Lf)(X t )], where Ex stands for the expectation with respect to P st 



conditioned on the initial state Xq, see e.g. [11, Ch. X]. The time reversed dynamics is 



generated by the adjoint L* of L with respect to the invariant measure p, that is 

ET\fLg] = E»[{L*f)g] (1.3) 

The measure p, which is the same for both processes, is a distribution over the configu- 
rations of the system and formally satisfies pL = 0. The expectation with respect to p is 
denoted by E^ and f,g are observables. We note that the probability P st , and therefore 
P* t , depends on the invariant measure p. The finite dimensional distributions of P s t are 
in fact given by 

Pst (X tl = • ■ • ,X tn = (f) tn ) = p((f) tl ) Pt 2 - tl (0 tl (f> t2 ) ■ ■ ■ Ptn-tn-A&tn-j -> <Pt n ) (1-4) 



where Pt{<pi —> <t>-z) is the transition probability. According to (|1.2|) the finite dimensional 
distributions of P* t are 

P st { X t, =&!>••• , X l = <Pt n ) = KMPti-hifal -*&!,)••• Pt-tn-lifan-l 0tn) 

= KfaJPU-U-li&tn -> <Pt n ^) ■ ■ ■ Pt 2 -tMt2 -» 0*1 ) (1-5) 

in particular the transition probabilities pt{4>\ 4>2) anc i Pti^i ~ * 4>2) are related by 

M0i)pt(0i -» 02) = Kfa)Pt{fc -> 00 (i- 6 ) 
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which reduces to the well known detailed balance condition if pt{4>\ — ► 4>2) = Pti&i ~> 4>2)- 
We require that also the evolution generated by L* admits a hydrodynamic description, 
that we call the adjoint hydrodynamics, which, however, is not necessarily of the same 
form as In fact we shall discuss a model in which the adjoint hydrodynamics is non 

local in space. 

In order to avoid confusion we emphasize that what is usually called an equilibrium 
state, as distinguished from a SNS, corresponds to the special case L* = L, i.e. the detailed 
balance principle holds. In such a case P st is invariant under time reversal and the two 
hydrodynamics coincide. 

4) The stationary measure P s t admits a principle of large deviations describing the 
fluctuations of the thermodynamic variable appearing in the hydrodynamic equation. 
This means the following. The probability that in a macroscopic volume V containing N 
particles the evolution of the variable p^ deviates from the solution of the hydrodynamic 
equation and is close to some trajectory p(t), is exponentially small and of the form 

Pst (p N (X t ) ~ p(t),t e [h, t 2 \) « e -"[WO)+ J [«i,t a ](/>)] = e - NI lH,H\(P) (1.7) 

where J(p) is a functional which vanishes if p(t) is a solution of ( |1.1|) and S(p(ti)) is 
the entropy cost to produce the initial density profile p(ti). We adopt the convention for 
the entropy sign opposite to the usual one, so that it takes the minimum value in the 
equilibrium state. We also normalize it so that S(p) = 0. The functional J(p) represents 
the extra cost necessary to follow the trajectory p(t). Finally p N (X t ) ~ p(t) means 
closeness in some metric and ~ denotes logarithmic equivalence as N — > oo. This formula 
is a generalization of the Boltzmann-Einstein. We set the Boltzmann constant k — 1. 

This paper is divided in two parts. In the first one we present a general fluctuation 
theory of SNS based on the hypotheses formulated above assuming the knowledge of the 
functional J(p). The main results are outlined below. 

1. The Onsager-Machlup relationship has to be modified in the following way: "In a 
SNS the spontaneous emergence of a macroscopic fluctuation takes place most likely 
following a trajectory which is the time reversal of the relaxation path according to 
the adjoint hydrodynamics" . 

2. We show that the macroscopic entropy S(p) solves a Hamilton-Jacobi equation 
generalizing to a thermodynamic context known results for finite dimensional Lan- 
gevin equations |U| 0, (22]] . 

3. For a large class of systems we obtain a non equilibrium non linear fluctuation 
dissipation relationship which links the hydrodynamic evolutions of the system and 
of its time reversal to the thermodynamic force, that is the derivative of S(p). If 
S(p) is known this relationship determines the adjoint hydrodynamics. 

4. From the last two results we derive an H Theorem for S(p): it is decreasing along 
the solutions of both the hydrodynamics and the adjoint hydrodynamics. 
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In the second part we test the theory outlined above in two boundary driven stochastic 
models of interacting particle systems: the zero range and the simple exclusion processes. 
The main results are outlined below. 

1. For the boundary driven zero range process, as observed in ||, the invariant mea- 
sure is a product measure. It is therefore possible to write the functional S(p), 
which in this case is a local functional of p, in a closed form and to construct the 
microscopic time reversed process explicitly. We derive both the hydrodynamics 
and the adjoint hydrodynamics. We compute the functionals J(p) and J*(p) and 
verify the generalized Onsager-Machlup principle and the fluctuation dissipation 
relationship. 

2. For the boundary driven simple exclusion process the invariant measure has long 
range correlations and is not explicitly known. The hydrodynamics has been ob- 
tained in j9|, [L(J; we obtain the asymptotics of the probability of large deviations, 
that is we calculate J(p). In one space dimension, Derrida, Lebowitz, and Speer |J 
have recently shown that the action functional S can be expressed in terms of the 
solution of a non-linear ordinary differential equation. We show how this result can 
be recovered by our approach: the Hamilton-Jacobi equation for S(p) (which is a 
functional derivative equation) can be reduced to the solution of the ordinary dif- 
ferential equation obtained in ||. By using the fluctuation dissipation relationship 
we also find the adjoint hydrodynamics. Moreover, in any spatial dimension, we 
can deduce a lower bound for the macroscopic entropy in terms of the entropy of an 
equilibrium state. In the one dimensional case this bound has been independently 
obtained in |J. 

Part of the results presented here have been briefly reported in |]J. Rigorous mathe- 
matical treatment of the boundary driven simple exclusion process will be given in ||. 

We conclude with some remarks to clarify the differences between equilibrium and non 
equilibrium states. The main problem in the SNS derives from the following situation. In 
equilibrium states the thermodynamic properties are determined by the Gibbs distribution 
which is specified by the Hamiltonian without solving a dynamical problem. On the 
contrary, in a SNS we cannot construct the appropriate ensemble without calculating 
first the invariant measure. At thermodynamic level, we do not need all the information 
carried by the invariant measure, but only the functional S(p) appearing in the generalized 
Boltzmann-Einstein formula ( |1.7| ). In general S(p), contrary to the equilibrium case, is a 
non local functional of the profile p. It turns out that the entropy S(p) can be obtained, 
both in equilibrium and non equilibrium, from J(p), which is therefore the basic object 
of the macroscopic theory. This step is simple for equilibrium, but highly non trivial in 
non equilibrium. 
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2. General theory 



2.1. Generalized Onsager-Machlup relationship 

We now derive a first consequence of our assumptions, that is the relationship between 
the action functionals I and I* associated to the dynamics L and L*. From equation 



L2[ ) and our assumptions it follows that P* t also admits a large deviation principle with 



functional I* given by 

with obvious notations. More explicitly this equation reads 

S(p(h)) + J( tlM (p) = S(p(t 2 )) + J { -t 2 ,- tl] (6p) (2.2) 

where p(ti), p(t 2 ) are the initial and final points of the trajectory and S(p(ti)) the entropies 
associated with the creation of the fluctuations p(tj) starting from the SNS. The functional 
J* vanishes on the solutions of the adjoint hydrodynamics. From ( j2.2|) we can obtain the 
generalization of Onsager-Machlup relationship for SNS. 

The physical situation we are considering is the following. The system is macroscop- 
ically in the stationary state p at t = — oo but at t = we find it in the state p. We 
want to determine the most probable trajectory followed in the spontaneous creation of 



this fluctuation. According to (177) this trajectory is the one that minimizes J among all 
trajectories p(t) connecting p to p in the time interval [— oo, 0]. From ( |2.2j ), recalling that 
S(p) = 0, we have that 

J [ -oo,o ] (p) = S(p) + J* , oo] (Op) (2.3) 

The right hand side is minimal if J* 0oo ^{Op) = that is if 9p is a solution of the adjoint 
hydrodynamics. The existence of such a relaxation solution is due to the fact that the 
stationary solution p is attractive also for the adjoint hydrodynamics. We have therefore 
the following generalization of Onsager-Machlup to SNS 

"In a SNS the spontaneous emergence of a macroscopic fluctuation takes place most 
likely following a trajectory which is the time reversal of the relaxation path according to 
the adjoint hydrodynamics" 

We note that the reversibility of the microscopic process X t , which we call microscopic 
reversibility, is not needed in order to deduce the classical Onsager-Machlup principle 
(i.e. that the trajectory which creates the fluctuation is the time reversal of the relaxation 
trajectory). In fact the classical Onsager-Machlup principle holds if and only if the 
hydrodynamics coincides with the adjoint hydrodynamics, which we call macroscopic 
reversibility. Indeed, it is possible to construct microscopic non reversible models in 



which the classical Onsager-Machlup principle holds, see O, [T3, [15 



2.2. The Hamilton- J acobi equation for the entropy 

We assume that the functional J has a density (which plays the role of a Lagrangian), i.e. 

J[t lM (p) = I 2 dt£(p(t),d t p(t)) (2-4) 
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From ( |2.3| ) we have that the entropy is related to J by 



S(p) = inf J[-ao,o]{p) (2.5) 



p 



where the minimum is taken over all trajectories p(t) connecting p to p. Therefore S 
must satisfy the Hamilton- Jacobi equation associated to the action functional J. Let us 
introduce the Hamiltonian TC(p,H) as the Legendre transform of C(p,d t p), i.e. 

H(p,H)= sup {(Z,H)-C(p,£)} (2.6) 

where (•, •) denotes integration with respect to the macroscopic space coordinates u, that 
is the inner product in L2(du). This notation will be employed throughout the whole 
paper. 

Noting that 7i(p,0) = 0, the Hamilton- Jacobi equation associated to ( |2.5|) is 

H(*£)=0 (2.7) 

This is an equation for the functional derivative A(p) = 5S/5p but not all the solutions 
of the equation TL (p, A(p)) = are the derivative of some functional. Of course only those 
which are the derivative of some functional are relevant for us. Furthermore, as well known 



in mechanics, the Hamilton-Jacobi equation (|2.7| ) has many solutions and we shall discuss 
later the criterion to select the correct one. 

Let also introduce the pressure G = G(h), where h = h(u) can be interpreted as a 
chemical potential profile, as the Legendre transform of the entropy S(p), namely 

G(h) = sup {(h,p)-S(p)} (2.8) 
p 

Then, by Legendre duality, we have 5G/5h = p and 5S/5p = h so that G{h) satisfies the 
dual Hamilton-Jacobi equation 

=0 (2.9) 



We now specify the Hamilton-Jacobi equation (|2.7|) for boundary driven lattice gases. 
We assume that the large deviation functional J may be expressed as a quadratic func- 
tional of d t p 

J[HM\{P) = \ f 'dt (V- I (^p-S(p)),x(p)- 1 V- 1 (^p-S)(p))) (2.10) 
= \ dt (v- 1 {dtp - \v ■ (D(p)Vp) ) , x(p)- 1 V- 1 [d t p - Iv ■ (D(p)Vp) ) ^ 

where D(p) = Dij(p) is the diffusion matrix in the hydrodynamic equation (|1 . 1|) and 
by V -1 / we mean a vector field whose divergence equals /. The form ( |2.10| ), which we 
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derive in the models discussed below, is expected to be very general; the functional J(p) 
measures how much p differs from a solution of the hydrodynamics and the matrix 
x{p) = Xi,j(p) reflects the intensity of the fluctuations. See II. 3.7] for a heuristic 
derivation of (|2.10| ) for reversible lattice gases. This form of J is also typical for diffusion 



processes described by finite dimensional Langevin equations | I2| . 

In this case the Lagrangian L is quadratic in d t p(t) and the associated Hamiltonian is 
given by 

H(p,H) = ~ (VH, x (p)VH) + ~ (H, V • (D(p)Vp)) (2.11) 
so that the Hamilton- Jacobi equation ( |2.7| ) takes the form 



We remark that the macroscopic entropy 5", given by the variational principle ( |2.5|) , 
depends only on the action functional J and is therefore stable with respect to microscopic 
perturbations which do not affect the dynamical large deviations. 

2. 3. The adjoint hydrodynamics and a non equilibrium fluctuation dissipation relation 

By assuming the quadratic form fl2.10|) also for J*, we deduce a twofold generalization 
of the celebrated fluctuation dissipation relationship: it is valid in non equilibrium states 
and in non linear regimes. 

Such a relationship will hold provided the rate function J* of the time reversed process 
is of the form ( |2.1CI|) with a different hydrodynamic equation (the adjoint hydrodynamics) 
that we write in general as 

d t p* = £*(p*) (2.13) 

with the same boundary conditions as 
We then assume J* has the form 



J(t lM (p) = \ f 'dt ((V- 1 (d t p - £*(/))) , xGT'V- 1 (d t p - ©•(/>))) (2.14) 
By taking the variation of the equation fl2.2|) , we get 



S(p)+D*(p) = V-(x(p)V^) (2.15) 

This relation holds for the non-equilibrium zero range process which we discuss later. We 
also note that it holds for the equilibrium reversible models for which the large devia- 



tion principle has been rigorously proven such as the simple exclusion process [19], the 



Landau-Ginzburg model and its non-gradient version |25|. It is also easy to check that 



the linearization of ( 2.15Q around the stationary profile p yields a fluctuation dissipation 
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relationship which reduces to the usual one in equilibrium. Accordingly, the matrix x{p) 
coincides with the Onsager matrix as defined in [|13|, [14], |27j . 

In order to verify the fluctuation dissipation relation ( [2.15 ), we need D(p), 2?* (p) and 
5S/5p. On the other hand, it can be used to obtain the adjoint hydrodynamics from 
D(p) and SS/Sp; the first is typically known and the second can be calculated from the 
Hamilton-Jacobi equation ( |2.12| ). 

Suppose we can decompose the hydrodynamics as the sum of a gradient of some 
functional V and a vector field A orthogonal to it in the metric induced by the operator 
K^ 1 where Kf = — V • (x(p)V/), namely 

with 

If 5V/Sp, like the thermodynamic force 5S/5p, vanishes at the boundary, it is easy to 
check that the functional V solves the Hamilton-Jacobi equation. 

Conversely, given S(p), by using the fluctuation dissipation relationship ( |2.15| ), we can 
introduce a vector field A(p) such that 



®(P) = \ V • (x(p) V^) + A(p) 



Sp, 

and the Hamilton-Jacobi equation implies the orthogonality condition 

Note the analogy with [|T^, Thm IV.3.1] for diffusion processes. 
2.4- H Theorem 

We show that the functional 5" is decreasing along the solutions of both the hydrodynamic 
equation ( |1.1| ) and the adjoint hydrodynamics 

d tP = ar(p) = v • (x(p)v^) - ©(p) (2.16) 

Let pit) be a solution of ( |1.1| ) or ( ^.16[ ), by using the Hamilton-Jacobi equation ( |2.12| ), 

we get 



piPit)) = (^(Pmd tP (t)) = -i(v^(p(t)), x(p(t))V^(p(*))} < (2.17) 
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In particular we have that ^S(p(t)) = if and only if |^(p(£)) = 0. Since we assumed 
there exists a unique stationary profile p, this implies that p is globally attractive also for 
the adjoint hydrodynamics ( |2.16 ). 

2. 5. A lower bound for the entropy S 

Let us consider any functional V(p), normalized so that V(p) = 0, whose derivative 
satisfies the Hamilton-Jacobi equation (|2.12|) and the condition 5V(p)/Sp = at the 
boundary. We shall prove the bound 

S(p)=mfJ hoofi] (p)>V(p) (2.18) 
p 

where the trajectory p(t) connects p to p. 

Fix ti < t 2 , two profiles pi, p 2 and a path p{t) in the time interval [ii, t 2 ] that joins pi 
to p 2 : p(ti) = pi, p{t 2 ) = p 2 - Let H, vanishing at the boundary, be given by the equation 

d tP = l -V ■ (p(p)Vp) - V • (x(p)Vtf) (2.19) 

We then claim that 

J[t lM (p) = V(p 2 ) - V( Pl ) 

Since the last term above is positive the bound ( |2.18| ) follows from the above identity. 
To prove ( |2.20|) we note, recalling (|2.10|) , that, since H is the solution of (|2.19| ), 

J^P) = I £ dt (VH, X (p)VH) 

We add and subtract in this expression V{5V(p)/Sp} to obtain 

We leave the first term unchanged and we show that the sum of the second and third 
gives V(p 2 ) — V(pi). Since 5V(p)/5p vanishes at the boundary, we may integrate by parts 
the second term; we get 

By the Hamilton-Jacobi equation ( |2.12| ), the third term is equal to 
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Adding together the previous two expressions, we obtain that the sum of the last terms 
in (|2.21|) is equal to 

jf* dt {iv • (D(p)Vp) - V • (x(p)VH) } 

Since p is the solution of ( j2.19| ), this expression is equal to 

* * dt (pjjr^p) = v(p(t 2 ))-v(p(t 1 )) = v(p 2 )-v( Pl ) 

which proves the claim. 

2. 6. Identification of the entropy 

In order to have a selection criterion for the solution V(p) of the Hamilton- Jacobi equa- 
tion, we consider the partial differential equation 

dtp = ~V ■ (D(p)Vp) + V • (x(p)V^) (2.22) 

As in the previous Subsection we assume V(p) is normalized so that V(p) = and satisfies 
SV(p)/8p = at the boundary. Note that this is the adjoint hydrodynamics Q2.16| ) if V 
coincides with S. 

If V = S we then have, by using the H Theorem fl2.17p , that the solution of the Cauchy 
problem (|2.22| ) with initial condition p relaxes to the stationary profile p so that 

\im V(p(t)) = V(p) =0 (2.23) 

Conversely if the above property holds, we can choose in (|2.20| ) the trajectory p(t) = p{—t), 
where pit) solves (f2~2l ). We then have H = 5V(p)/5p in ( ^2"0|) . The last term in ( ggtj ) 
becomes thus zero and V(pi) can be made arbitrary small; therefore ( |2.18| ) holds as an 
equality. 

The above argument shows that V — S if and only if ( |2.23| ) holds. 

2.7. Hamiltonian interlude 

As in Section 2.2, let us interpret the large deviation functional J as the action for the 
Lagrangian C(p,d t p), see ( |2.4|) , and J* as the action for the Lagrangian C*(p,d t p), see 
( j2.1-f ). Let also H(p,H) and H*(p,H) be the corresponding Hamiltonians obtained as 
Legendre transforms, see (|2.6| ). 

The time reversal relationship ( |2.2| ) implies the following relation between Lagrangians: 

C(p, dtp) = C*(p, -d t p) + (— , 8 t p) (2.24) 

As a consequence we obtain 

H(p,H)=H*(p, S -^-H^ (2.25) 
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Since C(p, d t p) and C*(p, —d t p) differ by a total time derivative, see (|2.24f) , we have the 
following. Given any p solution of either d t p = !D(p) or d t p = —D*(p) then p is a solution 
of the Euler-Lagrange equation for the Lagrangian C. Likewise given any p solution of 
either d t p = 2)*(p) or d t p = — 2)(p) then p is a solution of the Euler-Lagrange equation 
for the Lagrangian C*. 

In the case of the quadratic functional ( |2.10| ), we have 

£(P, dtp) = \ (V- 1 (dtp - D(p)) , X^-'V- 1 (d t p - D(p))> (2.26) 
The momentum conjugate to d t p is 

H = = -V- 1 ( X (P)- 1 V- 1 (d t p - D(p))) (2.27) 

where we recall that D(p) = |V- (D(p)Vp). Note that the above equation is the relation- 
ship Q2.19p ; as we shall see later, H is the external field we have to add to the microscopic 
dynamics to produce the fluctuation p. 

The Hamiltonian is given by ( 2.11 ), so that the Hamilton equations are 

9 tP = = 1 -V-(D(p)Vp)-V-( X (p)VH) 

dtH = = ~ [x[,(p)d Ut H d Uj H + D %3 (p)d Ui d U] H] ^ 

l<i,j<d 

where p(t, u) satisfies the same boundary conditions as in the hydrodynamical equation 
( P-.IP , H(t, u) vanishes at the boundary and x'(p) i s the derivative of x{p) with respect to 
P- 

We note that (p, 0) is an equilibrium solution of (|2.28| ) belonging to the zero energy 
manifold Tl(p, H) = 0. Any solution p(t) of the hydrodynamical equation ( |1 . 1| ) corre- 
sponds to a solution (p(t),0) of the Hamilton equation ( |2.28| ) which converges asymp- 
totically, as t — > +oo, to the equilibrium point (p, 0). The action J evaluated on this 
solution is identically zero; this corresponds to the trivial solution S = of the Hamilton- 
Jacobi equation (|2.12|) and is consistent with the vanishing of the conjugate moment 
H. Furthermore, if we take the time reversal of any solution of the adjoint hydrody- 
namics, i.e. dfp(t) = —D*(p(t)) we find a solution of the Hamilton equations given by 
(p(t),(5S/6p)(p(t))) which converges asymptotically, as t ^ — oo, to the equilibrium 
point (p, 0); the action J evaluated on this solution, as a function of the final state p, is 
the macroscopic entropy S(p). Both these trajectories live on the zero energy manifold. 
Similar properties hold for the Hamiltonian flow of TC*. 

Let us introduce the involution 9 on the phase space (p, H) defined by 

Q(p,H) = (p,^(p)-H) 

If we denote by resp. the Hamiltonian flow of Tl, resp. Tl*, then, by using (|2.25|) , 
as easy computation shows that acts as the time reversal in the sense that 

0$ t = $* t 9 (2.29) 
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equivalently, in terms of the Liouville operators, we have 



0{f,H} = -{Gf,H*} 

where / is a function on the phase space and {■, ■} is the Poisson bracket. 

The relationship (|2.29| ) is non trivial also for reversible processes, i.e. when TC = 7i*, in 
such a case it tells us how to change the momentum under time reversal. This definition 
of time reversal in a Hamiltonian context agrees with the one given in |23fl. 



3. Boundary driven zero range process 

We consider now the so called zero range process which models a nonlinear diffusion of 
a lattice gas, see e.g. [|T^]. The model is described by a positive integer variable r] x {r) 
representing the number of particles at site x and time r of a finite subset of the 
(i-dimensional lattice, A^r = Z d fl A^A where A is a bounded open subset of R d . The 
particles jump with rates g(r] x ) to one of the nearest-neighbor sites. The function g{k) is 
increasing and g(0) = 0. We assume that our system interacts with particle reservoirs at 
the boundary of A^v whose activity at site x is given by ip(x/N) for some given smooth 
strictly positive function ip{u). The microscopic dynamics is then defined by the generator 
(see for the one dimensional case) 

L N = L N)huik + L Nhound . 

where 

L NMlk f( V ) = \ E 9(v*)lf(v x ' y )-f(v)} 

(3.1) 



2 

x,y£A N 



Abound 



m = \ E {9(^)[f{v x n-f{v)]+^{y/N)[f{^ + )-f{v)}} 



\x-y\ = l 



in which 



r\ z if z^x,y 
■rffv = { r] z -\ if z = x (3.2) 
r] z + 1 if z = y 

is the configuration obtained from 7] when a particle jumps from x to y, and 

x,± _ f Vz if z ^x > 3 _ 

rj z ± 1 if z = x 

is the configuration where we added (resp. subtracted) one particle at x. Note that, since 
g(0) = 0, the number of particles cannot become negative. 

We also remark that if g (k) = k the dynamics introduced represents simply non 
interacting random walks on A^ (with the appropriate boundary conditions) in terms of 
the occupation numbers r) x . 
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3.1. Invariant measure 

Since the generator is irreducible (we can go with positive probability from any con- 
figuration to any other), under very general hypotheses on the function g(k) there exists 
a unique invariant measure. It is however remarkable that such invariant measure can be 
constructed explicitly (see |§ for the one dimensional case). 

Let Xn(x) be the solution of the discrete Laplace equation with boundary condition 
t/>, namely 

^A N X N (x) = i ^2 [^n(v) ~ X N (x)] = for any x G A N 

j/GZd (3.4) 
I x —3/|=l 

Xn(x) — ip(x/N) for any x ^ Ajy such that 3y G Ajy for which |x — y| = 1 

The invariant measure hn is the grand-canonical measure (an = ELeAjv obtained by 
taking the product of the marginal distributions 

Vx,n(Vx = k)= 1 V^) 

Z(Aiv(rr)) #(1) • • •#(&) 

where 

oo p. 

is the normalization constant. 

The fact that /i n is an invariant measure can be verified by showing that 

v 

for any bounded observable /. The above identity can be easily checked taking into 
account that, since Xjy solves ( |3.4j ), it is an harmonic function; in particular we have 

My) - A ^)] = ( 3 - 7 ) 

x£A N ,y£A N 
\x-y\ = l 

We emphasize that, if ip is not constant, the generator Ln is not self-adjoint with 
respect to the invariant measure so that the process is different from its time reversal and 
detailed balance does not hold. 

3.2. Hydrodynamic limit 

Let us introduce now the macroscopic time t = t/N 2 and space u = x/N. For u G A, 
t > 0, we introduce the empirical density as 



u ) = jjS E V*{N 2 t) S(u-^) (3.f 



13 



where 5 denotes the Dirac function. Note that 



/ dup N (t,u) = Yl Vx{N 2 t) 



is the average density of particles at (macroscopic) time t. 
Let G(u), u G A be a smooth function and consider 

(p N {t),G) = / dup N {t,u)G{u) 

J A 

To compute the time evolution of (pjsf(t),G) we first note that, according to the general 
theory of Markov processes, see e.g. fll], Ch. X], we have 

j f K (WW' G » = K {N 2 L N ( PN (t), G)) (3.9) 

where denotes the expectation with respect to the microscopic process with initial 
configuration rj. 

Let us assume that G has compact support K C A, so that only L^buik gives a non 
zero contribution; by summing by parts, we get 

N 2 L N (p N (t),G) = ~J^Y1 9MN 2 t))N 2 A N G(x/N) 

xgAm 

w \W £ 9MN 2 t))AG(x/N) (3.10) 

x£A N 

where we recall that A^v denotes the discrete Laplacian. 

At this point we face the main problem in the hydrodynamical limit: equation ( |3.10| ) is 
not closed in P7v(t) (its r.h.s. is not a function of pN(t)). In order to derive the macroscopic 
hydrodynamic equation we need to express g(r) x (N 2 t)) in terms of the empirical density 
PN{t)- This will be done by assuming a "local equilibrium" state (which can be rigorously 
justified in this context). Let us consider a microscopic site x which is at distance O(N) 
from the boundary and introduce a volume Q, centered at x, which is very large in 
microscopic units, but still infinitesimal at the macroscopic level. The time evolution 
in Q is essentially given only by the bulk dynamics L^buik; since the total number of 
particles in Q changes only via boundary effects and we are looking at what happened 
after 0(N 2 ) microscopic time units, the system in Q has relaxed to the canonical state 
corresponding to the density piy(t,x/N). 

Let us construct first the grand-canonical measure in Q with constant activity (p > 0, 
namely the product measure pg = YlxeQ with marginal given by 

1 if> k 

nth* = k ) 



Z(<p)g(l).-.g(k) 
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where Z{ip), which has been defined in ( |3.6| ), is the normalization constant. Let now Vq 
be the associated canonical measure at density a, i.e. 



^Vx = a\Q\ 
xeQ 



We introduce a function 0(a) by 



0(a) = \imE"%(g( Vx )) (3.11) 

Q-[Z d 

where we recall that E V Q denotes the expectation with respect to the probability Uq. 

According to the previous discussion, the system in the volume Q is well approximated 
by a canonical state with density pjv(i, x/N); we can thus replace, for N large, g(j] x {N 2 t)) 
on the r.h.s. of ( [3.1 OP by 0(ptv(£, x/iV)) thus obtaining 



±Eg{{p N (t),G)) « i<((0(p^(t)),AG)) (3.12) 



To see which are the boundary conditions satisfied by pjv(£) we need to consider the 
boundary dynamics, L^bound.- In contrast to the bulk dynamics, this is a non conservative, 
Glauber-like, dynamics. Since we are looking after 0(N 2 ) microscopic time units the 
density at the boundary has well thermalized to its equilibrium value which impose 

E^(<f)(p N (t,u)))^^(u) uedA (3.13) 

where the function has been defined above. 

Assume the initial configuration rj of the process is such that for any smooth function 
G we have 

lim (p N (0),G) = lim -L J" V*G(x/N) = [ d U1 {u)G{u) (3.14) 

xeAiv JA 

for some function 7. By the law of large numbers, pjv(t) becomes a deterministic function 
in the limit N — > 00, so that we can eliminate the expectation values in ( 3.12 ) and ( |3.13|) . 
In conclusion we have obtained, for any smooth function G, that 



lim (pjy(t),G) = / du p(t,u)G(u) 

N^oo ' J A 

where the convergence is in probability. Recalling ( |3.12| ), the limiting density p = p(t,u) 
solves 

d t p(t,u) = -Acj){p(t,u)) ueA 

<j)(p(t,u)) =ip(u) uedK ( 3 - 15 ) 
p(0,u) = 7(«) 
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which is the hydro dynamic equation for the boundary driven zero range process. 

We finally show that, by the equivalence of ensembles, we can express the function 
a i — > 0(a) introduced in ( 3.11 ), in a more convenient way, in terms of the grand-canonical 
measure \1q. By exploiting the product structure of pq and choosing the activity y?(a) so 

that pt^ a \r) x ) = a, we have 

0(a) = E# a) {g( Vx )) 

A straightforward computation shows now that <f)(a) = ip(at) so that the function a i— > 
0(a) is the inverse of the function ip \— > R((p) given by 



R(<p) = E*(r h ) = <p?M (3.16) 

i.e. R(<p) is the equilibrium density corresponding to the activity ip. From the assumptions 
on g it follows that R(<p) is strictly increasing. 

3.3. Dynamical large deviations 

In order to compute the probability of deviation from the typical behavior described by 
equation ( |3.15| ), namely the action functional J(p), we follow the classical strategy in large 
deviation theory: we need "only" to consider a perturbation of the system which makes 
the deviation p the typical behavior and write the probability in the unperturbed system 
in terms of the perturbed one. From this computation we shall extract, asymptotically in 
N, the factor exp{-N d J(p)}. 

We consider the zero range process in a (space time dependent) external field H(t,u) 
which is a smooth function of the macroscopic variables vanishing outside A, i.e. H(t, u) = 
for u ^ A, t > 0. The perturbed dynamics is specified by the time dependent generator 
Ln,t = -k/v,-r,buik + -^TV.r.bound. where 



x,j/eAjv 

+ ^y/N)e H ^ N2 ^- H ^ N2 ^[f(^ + ) - f(rj)}} 

The interpretation of the perturbed dynamics is the following: in the macroscopic scale, 
we simply introduced a small space-time dependent drift A^ _1 Vi?(t,M) in the motion of 
the particles. 

Assuming the initial configuration rj is associated to a density profile 7 in the sense of 
( p,14j ), by similar computations as the ones given for H = 0, we get that the hydrodynamic 
equation for the perturbed system is 

d t p(t, u) = ~A(f)(p(t, u)) - V • (0(p(t, u))VH{t, u)) u G A 

<j)(p(t,u)) = ip(u) uedK ( 3 - 18 ) 
p(0,u) = j(u) 
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If we regard pit, u) as given and consider H(t, u) as the unknown, the above equation 
tells us which is the perturbation for which p(t, u) is the typical behavior. Note (|3.18| ) is 
precisely the relationship ( 2 . 2 7] ) between dtp and the conjugate momentum H . 

Let us denote by P^, resp. f^ ,H , the probability for the unperturbed, resp. perturbed, 
process with initial configuration 77. We then have 



<^(A 2 t),t€[0,T]) 

dP»( v (NH),te[o,T)) = ex P{W^),#)} 

where 

J&MO.ff) = X) [HfrxMikitfT)- H(0,x/N)r) x (0)] 

r-T 

di ^ d t H(t,x/N)rj x (N 2 t) 

x£A N 

N 2 r T 



(3.19) 



2 



/ tft V g(v x (NH)) [ e H{t,y/N)-H{t, X /N) _ ^ 

Jo 



|as-y|=l 

A 2 ' T 



2 



/eft V ^/A) [ e ^(Wiv)-if(W^) _ 1] (3.20) 



x£A N ,y£A N 
\x-y\=l 



See Appendix |A] for a derivation of the above formula (another proof can be found in [ TE , 
Prop. Al.7.3]). 

With the help of equation ( |3.19| ) we can write the probability that pjv(^) is close to 
p(t) for the unperturbed system as follows 

K (PN(t) ~ * e [0, T}) = E»> H {e- J io^ H h PN ^ (3.21) 

where E^' H denotes the expectation with respect to the perturbed process. Using the 
explicit expression ( |3.2(J| ), by Taylor expansion, we get 

j£ T] ( V (-),H) » N d {(p N (T),H(T))-(p N (0),H(0))- F dt {p N {t),d t H{t)) 



\ [ dt Jji E ^M* 2 *)) [*H{t,x/N) + |V#(t,*/A)| 2 ] 



2 



xeAjv 

T 



<14 



+- / eft / dui){u)dnH(t,u)\ (3.22) 



where daH(t,u) is the normal derivative of H(t,u) (h being the outward normal to A). 

If r](N 2 t) is a typical trajectory for the perturbed process, by the same argument 
given in derivation of the hydrodynamical equation, we can replace g(rj x (N 2 t)) above by 
4>(pN(t,x/N)). Since pjv(t) is close to pit), and 

lim (p N (t) ~ p(t), t€[0,T\) = l 

N— >oo ' 
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from ( |3.21| ) we get 



(p N (t) ~ p(t), t e [0,T]) « exp {-iV d J [0iT] (p)} 



(3.23) 



where 



J M (p) = <p(T),fr(T)>-(p(0),fT(0)> 



eft (p(t),d t H(t)) 



-\l dt (<f)(p(t)),AH(t) + \VH(t)\ 2 ) + l- [ dt [ du^{u)d h H{t,u) (3.24) 

^ JO ^ JO J.9A 

Recalling that H and p are related by ( |3.18| ), we finally get, after an integration by parts 
(recall that H(t,u) vanishes at the boundary dA) 



J[o,T](p) = If dt <0(p(t)),|Vtf(t)| 2 ) 



(3.25) 



The action functional J is defined to be infinite if p does not satisfy the boundary con- 
ditions in ( |3.18| ). From ( |3.25| ) and ( |3.18| ) we see that J[o,t] is of the form ( |2.10| ) with 

A j = fl(j>)8ij and Xi,j(p) = Hp)h,j- 

The rigorous derivation of the action functional J requires some difficult estimates. 
In fact, while in the proof of the hydrodynamic limit it is enough to show that we can 
replace g(rj x (N 2 t)) by <p[p^(t,x/N)) with an error vanishing as iV —>■ oo, in the proof 
of the large deviations we need such an error to be o(e~ Nd ). This can be proven by the 
so called super exponential estimate, see [IS, ITPI, which is the key point in the rigorous 
approach. 

S.Jf.. Macroscopic entropy and adjoint hydrodynamics 

From the expression ( |3.25| ) for J it follows that the Hamilton-Jacobi equation ( |2.12| ) for 
the boundary driven zero range process is 



Sp 



Let us consider the functional 



S(p) 



du 



p{u) log 



A(«) 



log 



Z(<P(p(u))) 
Z(X(u)) J 



(3.26) 



(3.27) 



where Z{tp) has been defined in ( |3.6| ) and X(u) is the stationary activity profile, namely 
X(u) = (p(p{u)) where p is the stationary solution of the hydrodynamic equation (|3.15|) . 
Note that A is also the macroscopic limit of A at, solution of (|3.4p. By using that <p i— > R((p) 
given in ( p,16|) is the inverse function of p i— > <f)(p), we get 



SS(p) 
5p{u) 



\og(j){p{u)) - log A(m) 



(3.28) 
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We remark that the functional S given in ( |3.27| ) is uniquely characterized by fl3.28| ) once 
we impose the normalization S(p) = 0. 

An easy computation shows that the functional S given in (|3.27| ) solves the Hamilton- 
Jacobi equation ( |3.26|) . Recalling that (f)(p{u)) = \{u) = ip(u) for u G dA we have indeed 



V log M, mV log M\ + /log M, A0(,)\ 



V log ™ 
, W ,^_(v Wp )-A] ™)-( V A™ 

«H/>), ^\ - (VA, ^ + (W) - A] , V — ') = (3.29) 

since AA(tt) = for u G A. 

From the fluctuation-dissipation relationship (|2.15| ) we get that the adjoint hydrody- 



namic equation for the boundary driven zero range process is 

d t p*{t,u) = ~A4>(p*(t,u)) - V • (0(p*(*,u))VbgA(u)) u G A 

(j)(p*(t,u))=i>(u) uedA ( 3 - 30 ) 
p*(0,tt)= 7 («) 

Recalling that A(w) = (f>(p(u)), the density profile p is also a stationary solution of ( |3.30| ). 
Since > 0> the right hand side of (|3.30|) is dissipative; therefore we have that 

p*(t) — > p as t — > oo; so that we meet the criterion (|2.23|) . 

Since in this model we know explicitly the invariant measure pn we can check whether 
the predictions (|3.27|) on the macroscopic entropy and ( |3.30|) on the adjoint hydrodynamics 
of the general theory are correct. 

Given a smooth function h{u) let us introduce the pressure G{h) corresponding to the 
chemical potential profile h as 

G(h) = lim^logE^ (exp{N d (p N ,h)}) 

= i^LW* E HZ (X N (x)e h ^)- log Z(X N (x))] 

= I du [logZ (X{u)e h{u) ) - log Z(X(u))] (3.31) 

J A 



where Z(cp) has been defined in (|3.6|) and Xn(x) is the solution of (|3.4|). 
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By a standard computation due to Cramer we have that the Legendre transform of 
the pressure G(h), i.e. the macroscopic entropy, is the rate functional for the asymptotic 
probability of large deviations of the density profile in the invariant measure p^. Let us 
in fact introduce a perturbed measure p h N by 



»n(v)= II Z (\ N (x)e h (*/ N )) expfoOE/W)*?*} fJ-xAVx) 



_Z{\ N {x)) 

xgAjv 

and, for a fixed density profile p(u), choose the chemical potential profile h so that 

E^(rj x )=p(x/N) (3.32) 
namely h(x/N) = \og[4>(p(x/N)) / X N (x)}. We then have 

MPN-p) = ( e -E, 6Aw [M^)^-logX(A J v(x)eMViV) )+logZ(Ajv(x))]l ^^^ 

« fe-^^'^-^l^^) « e -^lO».fc>-GW] (3.33) 



since, by the law of large numbers, lim jv->oo P%(pn ~ p) — 1- From ( |3.32| ) we have that 



5G/5h = p. We therefore have obtained precisely that 5 is the Legendre transform of G. 
A computation, which is left to the reader, shows now that the Legendre transform of the 
right hand side of ([3.31|) gives indeed the functional S(p) defined in (|3.27| ). 



We want to stress a main difference between the macroscopic computation fl3.29|) and 



the microscopic one just given. While the former depends on the action functional J, 
which involves only macroscopic quantities, the latter depends on the explicit expression 
of the invariant measure p^. In particular the macroscopic computation is independent 
of the specific way the interaction with reservoirs is modeled (provided of course the 
functional J is not affected). 

We now discuss the adjoint hydrodynamics from a microscopic point of view. Since 
the invariant measure p^ is explicitly known we can obtain the adjoint generator L* N 
which is defined by the identity (|L3|). Recalling that solves (|3.4j) and that ( p. 71) holds, 
we have that L* N = L* N hulk + L* Nhound where 

WW = I E 9(v*)yH[f(v x ' v )-m] (3-34) 

\">— 1/1=1 

^bound./W = \ E {^(^^^[/(^~)-/(^)]+A J v(x)[/(^ + )-/(^)]} 



|as— vl=l 



Notice that the form of (|3.34|) is the same as ( p.l| ) with the rates modified in such a way to 



invert the particle flux. The generator L* N can also be interpreted as the original system 
perturbed by a time independent external field H(t,x/N) = log Ajy(x), compare ( |3.17| ) to 



( [3.34! ). In particular we have that the adjoint hydrodynamic equation is indeed given by 
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By repeating the same argument given in Section 3.3, it is not difficult to show that the 
action functional J* describing the probability of large deviations from the hydrodynamic 
behavior for the adjoint process is given by 

J[0,T](P) = \[ dt WW)' l Vi ^)i 2 ) ( 3 - 35 ) 

where K(t) = K(t, u) is to be obtained from p by solving the equation 

f d tP (t, u) = -A<f>(p(t, u)) - V • (4>(p(t, u))V [log X(u) + K(t, u)}) ueA 

\ (j>(p(t,u))=il>(u) uedA ( 3 - 36 ) 

A computation now allows us to check that the identity ( |2.2|) , which has been obtained 
only by a time symmetry argument, holds for the boundary driven zero range process. 



4. Boundary driven simple exclusion process 

The simple exclusion process is a model of a lattice gas with an exclusion principle: a 
particle can move to a neighboring site, with rate 1/2 for each side, only if this is empty. 
Let, as in the previous Section, Ajv = 1> d H iVA and denote by T] x (t) e {0, 1} the number 
of particles at the site x at (microscopic) time r. The system is in contact with particle 
reservoirs at the boundary of A^ whose activity at site x is given by ip(x/N) for some 
given strictly positive smooth function i[)(u). 

The microscopic dynamics is defined by the generator = Ljv.buik + -Abound, where 

LNWkfiv) = \ E V*(l-Vy)[f(v x ' y )-m] (4-1) 

\x-y\=l 



J N,bonnd 



-fW = \ E {v4f(v x n-m]+^(v/N)(i-v x )[f(v x ' + )-f(v)]} 



\x-y\=l 



where r] x,y and r] x,± have been defined in (|3.2j) and (|3.3|) . 

In contrast to the zero range model the invariant measure pj^ is not known explicitly; 
we shall see that it carries long range correlations making the entropy non local. 

4-1. Hydrodynamic equation and dynamical large deviations 

The hydrodynamic equation for the simple exclusion process can be derived by the same 
argument given for the zero range process; in fact it is easier in this case because a com- 
putation analogous to (|3.10|) leads directly to a closed equation in the empirical density. 
We find that the limiting density evolves according to the linear heat equation 

d t p(t, u) = -Ap(t, u) u e A 
p(0,u) = j(u) 
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where 7 is the initial density profile, associated to the initial microscopic configuration 
f] in the sense ( [3.14|) . In this case the density of particles p takes value in [0, 1]. The 
hydrodynamic limit for more general boundary driven models has been discussed in P, 
As in the previous Section we shall denote by p = p(u) the unique stationary solution of 
C2D. 



The action functional J describing the probability of large deviations from the hydro- 
dynamic behavior can be obtained as for the zero range process. In this case the perturbed 
dynamics is defined by the time dependent generator Ljy = L^ T bulk + -^^ iTi bound. where 



x,y£A N 
\x-y\=l 

L H N , T , M .m = \ y, {^f**™-™**™ [ftfn - m] ( 4 - 3 ) 



and the external field H(t, u) vanishes for u ^ A. 

The hydrodynamic equation for the perturbed dynamics is then given by 



d t p{t, u) = -Ap(t, u) - V • |/3(t, u) [1 - p{t, u)\ VH(t, u) } u e A 

**•»> = IT,' 

k p(0,u) = j(u) 



1 + ip(u) 



For the simple exclusion process the functional J^ T ^ defined in ( |3.19| ) is given by 
J& T] (ri(-),H)= [H(T,x/N) Vx (N 2 T) - H(0,x/N) Vx (0)] 

dt d t H(t,x/N)r) x (NH) 
v- r T 

~ dt ]T Vx(N 2 t)[l-rj y (N 2 t)] [ e H(t,y/N)-H(t, x /N) _-q 

J ° x,yeA N 
\x-y\=l 

~ fdt Yl VxiNH) [ e mt,y/N)-H(t,x/N) _ X ] 



x<=A N ,ygA N 
\x-y\=l 

t2 rT 



M 2 r 1 

~ dt ^V/ N ) i 1 - V*(N 2 t)] [e H ^- H ^ - 1] (4.5) 

J ° x&A N ,y£A N 
\x-y\=l 

we refer again to Appendix |A] for the derivation of the above formula. 
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By Taylor expansion, summation by parts, the replacements f] x (N 2 t)[l — r] y (N 2 t)] rs 
p N (t,x/N)[l - p N (t,x/N)] in the bulk and r] x (N 2 t) « ip(t,x/N)/[l + ip\t,x/N)\ at the 
boundary (which can be heuristically justified by the same argument given for the zero 
range process) we get 



J^ T] ( V (-),H) » N d [(p N (T),H(T)) - (p N (0),H(0)) 

dt L N (t),d t H(t) + ~AH(t)\-±£dt (p N (t){l-PN(t)},\VH(t)\ 2 ) 



+1 f dt [ du - , dnH(t,u) } (4.6) 



2 Jo J ah 1 + 

Recalling that the hydrodynamic equation for the perturbed dynamics is given by 
after an integration by parts, we finally get the action functional J 



1 



T 



J[o,T](p) = 2 I dt (P^ 1 " > l VJ W) ( 4 - 7 ) 







which is of the form ( |2.10|) with D it j(p) = 5 it j and Xi,j(p) = p[l — p]Sij- As for the zero 
range, the functional J[o,r](p) is defined to be infinite if p does not satisfy the boundary 
condition in ((4.41). 



Reduction of Hamilton- J acobi to a non linear differential equation (d = 1) 
We consider here the boundary driven simple exclusion process in one space dimension. 
In a very interesting paper, by using a matrix representation of the microscopic invari- 
ant measure and combinatorial techniques, Derrida, Lebowitz and Speer have recently 
shown that the action functional S (which we called the macroscopic entropy) can be 
expressed in terms of the solution of a non-linear ordinary differential equation. We 
show next how this result can be recovered in our approach by following the dynami- 
cal/variational route explained in Section Namely, we consider the variational prob- 
lem Q2.5p for the one-dimensional simple exclusion process and show that the associated 
Hamilton-Jacobi equation which, taking into account (|4.7|) and ( |4.4|) , is the functional 
derivative equation 

V— ,p(l -p)v— j + (j^> a p) = ( 4 - 8 ) 

can be reduced to the non-linear ordinary differential equation first obtained in ||. 

For notation simplicity, we assume that A = (—1,1), so that dA = ±1. We shall 
also assume the macroscopic density profile p = p(u) satisfies the boundary conditions in 
equation ( |4.2j ). 

We look for a solution of the Hamilton-Jacobi equation (|4.8|) by performing the change 
of variable 

log , n , -<P(u;p) (4.9) 



Sp(u) 1 - p(u) 
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for some functional (f)(u; p) to be determined satisfying the boundary conditions 0(±1) = 
logp(±l)/[l-p(±l)]=logV(±l). 

Inserting (|4.9|) into (|4.8|), we get (note that p — / (1 + e^) vanishes at the boundary) 



= -/v(log^-0j,p(l-p)V0 
= -(Vp,V0) + <p(l-p),(V0) 2 > 



We obtain a solution of the Hamilton-Jacobi if we solve the following ordinary differential 
equation which relates the functional <p(u) = <f>(u\ p) to p 



-p(«) we (-1,1) 



[V0(m)] 2 1 + e^) 

0(±l)=log^(±l) 
A computation shows that the derivative of the functional 



(4.11) 



S(p) 



j du jp log p + (1 - p) log(l - p) + (1 - p)0 - log (1 + e*) + log ||| (4.12) 



is given by ( OP when 0(w; p) solves ( |4.11|) . According to the discussion in Section |2~6| , 
we will be able to conclude that Q4.12| ) is indeed the macroscopic entropy as soon as we 
show that it meets the criterion (|2.23|) . This will be done in the next Subsection. By the 
change of variable = log[F/(l — F)] equation ( [4. lip becomes the one obtained in 0. 

One may be tempted to repeat the same computation in arbitrary dimension; one 
would obtain a partial differential equation analogous to ( 4.11Q . However, in more than 
one dimension it does not exist, in general, a functional S whose derivative is given by 
( [4.9D with cf) and p related by such partial differential equation. 

The equation ( |4.11| ), considered as a relationship expressing p in terms of 0, can be 
interpreted in the following way. Let 

S (p) = J <Mplogp+(l-p)log(l-p)} (4.13) 

be the equilibrium entropy. Since 6So/5p = log[p/(l — p)] we have 

5S 5S 

nu-p) = - — — 

op op 

If G(<f>) is the Legendre transform of So — S, we find that 5G/5<ft = p which is the rela- 
tionship (|4TTl| ). 
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We note that the remark after ( |3.33| ) for the zero range process also applies to the 
present context. In particular the macroscopic computation ( [4.1 0|) depends only on the 
action functional J and is therefore stable with respect to microscopic perturbations which 
do not affect the dynamical large deviations. 

4-3. Adjoint hydrodynamics (d = 1 ) 

By using the fluctuation dissipation relationship (|2.15| ) and the expression (|4.9| ) for 5S/ 5p, 
we find that the adjoint hydrodynamics is given by the equation non local in space 



d t p*(t,u) 

^' ±1): 1 + *( ± 1) 
[ p*(0,w)= 7 («) 



-Ap*(t, u) - V { P *(t, u)[l - p*{t, u)]V(j>{u- P \t))} u e (-1, 1) 
>(±1) 



(4.14) 



where 0(w; p*(t)) is to be obtained from p*(t) by solving ft4.11|) . Since p) = log[p/ (1 — 
p)], we see that p is also a stationary solution of ( |4.14j ). 

We now show how (|4.14j ) can be related to the heat equation. Let p*(t,u) be the 
solution of ( |4.14| ) and introduce F = F(t,u) as 



F(t,u) 



(u;p*(t)) 



(4.15) 



1 _|_ e <t>(u;p*(t)) 

it is not too difficult to check (see Appendix |FJ) that F(t, u) satisfies the heat equation 

1 



d t F(t,u) 
F(t,±l) 

F(0,u)-- 



■AF(t,u) ue (-1,1) 
^(±1) 



i + v(±i; 



(4.16) 



X _|_ e <t>(u-n) 

Conversely, given F = F(t,u) which solves ( [4.16D , by setting 

AF(t,u) 



p*(t,u) = F(t,u) + F(t,u)[l - F(t,u)) 



[VF(t,u)f 



(4.17) 



a computation (see again Appendix |B|) shows that p* solves ( |4.14|) . 

We have thus shown how a solution of the (non local, non linear) equation ( 4.14 ) can 
be obtained from the linear heat equation by performing the non local transformation 
( |4.15| ) on the initial datum. In particular, since the solution F(t, u) of (|4.16|) converges as 
t — > oo to p, we see that the functional S(p) given in (|4.12| ) satisfies the criterion ( |2.23| ) 
so that it is indeed the macroscopic entropy. 

4-4- Non perturbative lower bound on the macroscopic entropy (d > 1) 

We discuss here a non perturbative bound for the macroscopic entropy in the bound- 
ary driven simple exclusion process in arbitrary space dimension d. Let Sq(p) be the 
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equilibrium entropy as defined in (|4.13 ), we shall obtain the following lower bound on 
S(p) 



S(p) > S (p) - S Q (p) - (p- P.^(p) 
du ip log - + (1 — p) lor 



P 



l-p 



S(P) 



(4.18) 



with a strict inequality for p ^ p. In the one dimensional case the bound ( 4.18 ) has been 
independently obtained in ||. 

Recalling that the dynamical action functional J of this model is given by fl4.7|) , a 
somewhat lengthy but straightforward computation gives 



J[-T,a]{p{-)) 



\ f_ T dt < v " - H . i dtp - \** 



S(p(0)) - S(p(-T)) 
+ 



( W ) + lAp - V • - p]V log ^) ) 



1 

+ 2 



-t A [p»(i-p-(«))] 2 



(p(t,u) -p(w))' 



(4.19) 



The last two terms on the right hand side of ( [4.19[ ) are positive. Therefore, if pit) is 
trajectory connecting p to p, we have 

S(p) = inf J[_oo,o](p) >3{p)- S{p) = S{p) 
p 

Moreover, since the last term on the right hand side of (|4.19 ) is strictly positive as soon 
as p 7^ p, we have the strict inequality in ( |4.18|) for p ^ p. 

4-5. Perturbative solution of the Hamilton-, J acobi equation (d>l) 

We show here how the Hamilton- Jacob i equation ( 2.12j ) can be used to get a perturbative 
expansion for the entropy S around the stationary profile p. We discuss only the expansion 
up to the second order but it will be clear how to generate an iterative approximation 
scheme. 

From a computational point of view it is convenient to expand the pressure G(h) 
defined in (|2.8|). Since p(u) = SG(h)/Sh(u) and h{u) = 6S(p)/6p(u), the dual Hamilton- 
Jacobi equation (|2.9|) in this model is 



V/i, 



SG 



Sh 



5G 



5h 



Vh 



„, ^SG 



(4.20) 



where h(u) satisfies the boundary conditions h(u) = for u G dA. 
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Recall that p(u) is the stationary solution of the hydro dynamic equation ( |4.2| ) and 
introduce 



G(h) = I du log [1 + p(u) (e h{u) - l)] 



(4.21) 



Note that G(h) is the Legendre transform of S(p) defined in (|4.18| ). We look for a solution 
of ( |4.20| ) in the form 



(4.22) 



G{h) = G(h) + (g, h) + -{h, Bh) + 0(h 3 

for some function g = g(u) and some linear operator B. 
Note that S(p) has a minimum for p = p and 



s(p) = -(p-p,c- 1 (p-p)) + o ((p - pf) 

where C is the covariance of the density fluctuations with respect to the invariant measure. 
Therefore 



G(h) = (p,h) + ^(h,Ch) + 0(h 3 ) 



(4.23) 



hence the second derivative of G at h = is the covariance of the density fluctuations. 
By comparing ( |4.22|) to ( |4.23|) we get 



C 



5 2 G 



Sh 2 



h=0 



B = p(l-p)l + B 



(4.24) 



Since G is the pressure for the equilibrium system the operator B represents the contribu- 
tion to the covariance due to the non equilibrium boundary conditions. For the boundary 
driven simple exclusion process the covariance of the fluctuation has been derived in [||. Q 
where it is shown that it is non local in space. Therefore the perturbative expansion of 
the Hamilton-Jacobi equation will give a different derivation of the covariance. 
We have 





sg\ 


Sh \ 




5G 
~Sh ~ 


P(l- 


[1+P( 



p(l - p)e 



[l+p{e h -i)X 
1 



[l + p(eh-l)} 



so that by plugging ([4. 22|) into ( |4.2(]| ) and expanding up to second order in h we get 

(V/i, [1 + (1 - 2p)h] Vp + V# + VBh) = 
Recalling that h vanishes at the boundary, we thus get g = and 

v(^Y(l-2p)Vp\ = (h,ABh) 
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which, after an integration by parts, becomes (recall that Ap = 0) 



(h,ABh) = (h,\Vp\ 2 h) 



The operator B therefore satisfies 



~ [AB + BA] = |Vpf 



In particular, if Vp is constant, B has the kernel 

B{u,v) = iVppA- 1 ^,^) 



(4.25) 



where A -1 (w,v) is the Green function of the Laplacian (with Dirichlet boundary condi- 
tions). The fact that B is a negative operator can also be obtained as a consequence of 
the bound S(p) > S(p). 

By analogous computation one can obtain also the higher order terms in the expansion 
of the pressure which are the higher order cumulants. In the one dimensional case Vp is 
constant and we state below the result of the expansion up to the third order. 



A. Derivation of flOOp and (fO|) 

Let Q be a finite set and consider a continuous time jump Markov process X t on the state 
space Q with generator given by 



where the rate A is a positive function on Q and p(u),u') is a transition probability. We 
can construct a realization of X t as follows. Fix an initial condition X = uj . The process 
waits an exponential time T\ with rate A(co> ) and then jumps to uj\ with probability 
p(cuo,uJi)] the law of T\ is 



The process waits now an exponential time T2, independent of n, with rate A(u;i) and 
then jumps to uj 2 with probability p{ijj\^uj2), and so on. Consider the piecewise constant 
trajectory X s , s G [0,T] with n jumps given by 





(A.l) 




(A.2) 



< s < ti 

h < s <t 1 +t 2 



X s 



= < 



(A.3) 



Un-l ti+t 2 -\ + < s < t x + t 2 H h t 

u n t\ + h H + t n < s < T 



■n 
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Its probability density is given by 

n 

dP„ (X s , s G [0,T]) = J] (A(o; i _ 1 )p(o; i _ 1 ,a; i ) e -*<<*-0*< dt^j ■ e -K^-^] 



i=i 



where a n = t\ + ■ ■ ■ + t n 

If A and p depend explicitly on time we can construct a realization X t by the same 
procedure; in such the law of T\ is 



p( Tl <t)= ds A (wo, s) e -Iods'X^ ,s') (A.4) 
and analogous distributions for Tj. We thus get 

n 

dP^ (X s , s G [0,T]) = J] (A^,^^^,,^;^) dU) . e -^dsHX s , s) (A 5) 



where = ii + • • • + tj. (resp. a = 0) are the jump times of X s . 

Let us now consider two processes X t (resp. X' t ) of this type with rates X(u,t) (resp. 
A'(u;, £)) and transition probability p(co, u/; t) (resp. p'(cu, u)'\ tj). We can write the formula 
( |A.5| ) also for the process X' and, by taking the ratio of those expressions, we get the so 
called Radon-Nikodym derivative 



dP' 

-fHX s , s G [0,T]) = exp {J [0> T](X)} 



= exp <^ V log <T - 1 ' 7 / / - / A' (X a , s) - \(X a , s 

\ti \{X ai _ l ,a i )p{X (Ti _ 1 ,X ai - 1 a i ) J L 1 ' 1 

We now consider the special case in which X has generator given by (|A.1|) and X' = X F 
has a time dependent generator 

Lff(co) = A(w)p(w, wO^*)-^*) [/(a/) - /(a;)] (A.7) 
a/ en 

which is of the same form with 

X'(co,t) = Ku)p{u,oo')e F{ -"'ti- F ^ 

TtMt) = -j i-AMM^Oe^-^ 
A '{u, t) 

From ( |A.6|) we get that the Radon-Nikodym derivative is given by 

dP F 

-p^(X s , s G [0,T]) = exp{J [0tTi (X,F)} 
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with 



J [0 , T] (X,F) = [F(X ai , o-j) — F(X t7i _ 1 ,a, 

i=l 

- [ T dsX(X s )J2p(X s ,oo') \fVA-nx.A _ ! 



'en 



£ 

i=l 



ds d s F(X ai _ 1 ,s) 



(A.8) 



- / dse~ F ^Le F ^ 
Jo 



= F(X T ,T)-F(X ,0)- I ds[d s F(X s ,s) + e- F{x ^Le F( - x ^} 

Jo 

Formulae (|3.20|) and ( fl5|) are special cases of (|A.8|) obtained by choosing 

F(tj,t) = H(t/N 2 ,x/N) Vx 



x&Ar 



B. Adjoint hydrodynamics for one dimensional simple exclusion 

Let p*{t) be a solution of the adjoint hydrodynamics for the one dimensional simple 
exclusion process ( f4.14|) . By the remarks in Subsection pT7| , (p(t),H(t)) with 

p(t) = p*(-t) 



(B.l) 



is a solution of the Hamilton equations ( |2.28p which for this model read 

dtp = l -Ap-v[p{l-p)VH) 
d t H = - l -{l-2p)(VH) 2 - l -AH 



(B.2) 



By plugging (|B.1| ) into ( |B.2| ) and performing the change of variable ( |4.15| ), a straightfor- 
ward computation yields 



dtp* 



d t F 



1 (VF) 2 
— -AF + (p* - F)—. 

2 KH J F(1-F) 



(B.3) 



By writing the equation ( }4.11|) in terms of F = e^/ (1 + e^) and p replaced by p* we get 



(p*-F)=F(l-F)- 



AF 



( Vi T 

F(t,±l) = p*(t,±l) - 



for any (t, u) G [0, oo) x (—1, 1) 



(B.4) 



i+<M±i) 
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which inserted in ( |B.3| ) concludes the proof that F(t, u) as defined in ( |4.15| ) satisfies the 
heat equation. 

The converse statement, namely that if we define p* = p*(t,u) as in fl4.17| ) (with 
F = F(t,u) the solution of (|4.16|) ) then it satisfies the non local equation (^4.14|) , can be 
checked without invoking the Hamiltonian formalism. Indeed, from ( [4.1 7| ) we get that 



^ (1 -^-l + (l-2F)^-Fa 



F(l-F) 



(VF)' 



(AFf 
(VFY 



(B.5) 



recalling (|4.16|) , by a somehow tedious computation of the partial derivatives which we 
omit, we get 



d t A 

2 



F(l-F) 



AF 

(vf7 



F(l-F) 



(B.6) 



Therefore, recalling ( (4.17|) , the function p*(t) satisfies 



d t p* = -Ap* - V 



p*(l-p* 



-VF 



(B.7) 



F(l-F) 

which is precisely ( f4.14| ) written in terms of the variable F = F(p*) instead of <fi = <fi(p 
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